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We introduce a method of quantum tomography for a continuous variable system in position 
and momentum space. We consider a single two-level probe interacting with a quantum harmonic 
oscillator by means of a class of Hamiltonians, linear in position and momentum variables, during 
a tunable time span. We study two cases: the reconstruction of the wavefunctions of pure states 
and the direct measurement of the density matrix of mixed states. We show that our method can 
be applied to several physical systems where high quantum control can be experimentally achieved. 
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I. INTRODUCTION 

The tomography of low-dimensional systems [T] is a 
well established process in fundamentals of quantum 
physics and quantum information [21. Measuring a fi- 
nite set of observables on identical copies of the system, 
one is capable of reconstructing the density matrix of a 
quantum state with a finite precision. However, practical 
experiments have a limited accuracy and the number of 
measurements grows exponentially with the size of the 
Hilbert space. Furthermore, such schemes do not have a 
direct translation to setups where the degrees of freedom 
are continuous, as is the case of atomic beams [3l|4], har- 
monic oscillators in trapped ions [5HTT] , cavity QED [T2l - 
[T3] , circuit QED P"5VE5] , nanomechanical resonators [21] , 
and superconducting qubits pS^BS] . 

In the case of continuous variable systems, the mea- 
surement of an infinite number of observables would be 
required and, in consequence, other simplifications or 
schemes have to be developed [211 [3D]. In linear op- 
tics or in atomic ensembles, where interactions are typ- 
ically quadratic, one assumes that the generated states 
are Gaussian. That is, they are fully characterized by 
first and second field quadrature moments, the so called 
covariance matrix |31| . In this case, again, a finite set of 
measurements suffices to determine the state, its entan- 
glement properties, and all other observables. 

However, not all physically realizable states are Gaus- 
sian. For example, photon-added states and superpo- 
sitions of Fock states [21] are interesting resources for 
quantum information processing. For a continuous vari- 
able description of such non-Gaussian states the most 
popular solution is the reconstruction of the Wigner func- 
tion [32 , a quasi-probability distribution that contains 
the same information as the density matrix. 

In this paper, we propose a method for reconstructing 
the state of a continuous variable system in position and 
momentum space, be for a pure state or a mixed state. 



To accomplish this purpose, we couple the quantum sys- 
tem to a two-level probe with an adequate interaction. 
For other reconstruction methods using different interac- 
tions and experimental requirements, see, for example, 
Ref. [33]. We show that, by monitoring the probe, we 
can obtain enough information to reconstruct the wave 
function of the system in position and momentum space. 
This paper is organized as follows. In Sec. II, 
we present the model Hamiltonian and we develop the 
method to obtain the wave function when the system is 
in a pure state and a mixed state. In Sec. Ill, we propose 
a physical implementation of our method. In Sec IV, we 
present our concluding remarks. 



II. QUANTUM TOMOGRAPHY OF 
CONTINUOUS- VARIABLE SYSTEMS 

From an operational point of view, our proposal re- 
quires the generic interaction Hamiltonian between a 
two-level probe and a ^-dimensional system 



H = hgaft ®(d-R + i3-P 



(1) 



Here, CTfj is a probe Pauli operator along n, a and /3 e M^ 
contain control parameters, while R — {Xi, . . . , Xg) and 
P = {Pi, . . . , Pi) are the dimensionless position and mo- 
mentum operators of the system. When the probe is cou- 
pled to the system, they exchange quantum information 
and, by monitoring the state of the probe as a function 
of time, we will be able to reconstruct the wave function 
i'siR) of the oscillator system. 

In the following we will consider the reconstruction of 
a quantum state, be it in a pure or a mixed state. For 
simplicity and without loss of generality, we will develop 
the theory for one-dimensional systems, £ = 1 in Eq. (Ill) 
and present useful examples. 



A. Pure states 

In many experimental setups suitable for quantum to- 
mography, the measurable quantity is the probe observ- 
able Gz] this could be for example a trapped and laser- 
manipulated two-level atom or a superconducting qubit. 
In this sense, the value of |'0s(a;)P will be obtained by 
monitoring the evolution of (tr^) [M]. To achieve this, 
we set the control parameters such that Eq. (fTl) turns 
into H — hgaUxX, where ax is the Pauli matrix along 
X direction and X = [a + a^)/-\/2 is the dimension- 
less position quadrature. The evolution of the two-level 
probe can be calculated as {(Jz)t = {'^\U'^<7zU\^), where 
U — exp{—iHt/h) and |^) is the total probe-system 
wavef unction. It is convenient to write 

{az)t - (*| cos ikX)az\^) + (*| sin {kX)ay\^), (2) 

with k ~ 2gat. From this expression, we are able to ob- 
tain |7/;s(a;)p by considering two different initial states for 
the probe, leading to two sets of measurements. In the 
first one, we consider the probe decoupled from the sys- 
tem such that we have an initial state |^) = jf/'s) (X" | t)2, 
such that, cTzl t)z = I t)z- In this case, the expectation 
value of az, Pz{k), will be 



P^ik) = f 

J —i 



COS {kx)\ip s{x)\^ dx. 



(3) 



We observe that the quantity Pz{k) is related to the even 
part of IV'sCa;)!^- Then, the second set of measurements 
provides us the corresponding odd part. This is found by 
preparing system initially in the state |\1/) = \ips) 'S) \ t)v 
such that CTyl t)i/ = I t)y In this case, the expectation 
value of az is given by 



P°{k) = / sm{kx)\il;s{x)\^ dx. 



(4) 



Combining the results of the two sets of measurements, 
we obtain the squared absolute value of the wave function 
by noticing that 



\A{x)\'^J'-'[P'zik)+^P:{k)] 



(5) 



where J^ ^[f{k)] = jw /-oo ■^(^)^ ^'^'^rffc is the inverse 
Fourier transform of /(/c) . In this manner, we have shown 
that the squared absolute value of the system wavefunc- 
tion can be obtained by monitoring the populations of 
a two-level probe. When the state of the system has no 
imaginary components, the measurement of |?/'s(a;)p is 
enough to obtain the wave function ips{x). This method 
for obtaining the wave function can be extended to the 
case of multimode quantum systems. 

To reconstruct the complete wavefunction when the 
state of the system has complex components, we will need 
additional information. In this case, we must properly set 
the control parameters to produce the Hamiltonian 



H = hgax{aX + (5P), 



(6) 



with tunable a and /?, and where P = i{a'' — a)/y/2 is 
the dimensionless momentum quadrature. Similar to the 
previous case, we can obtain 



P^^{4's\cos{kX + qP)\,p,), 
P: = {^bs\sm{kX + qP)\^s), 



(7) 
(8) 



where k — 2gta and q = 2gtf3. Using the Baker- 
Campbell-Haussdorf formulas, we can write 



P:^l{^Ps\e'''''+''^''+il.c.\^Ps 
P°^lm{^s\e'''^e"'Pe'''^^^\t^s] 



(9) 



(10) 



We consider now (qP) <C 1 and expand the previous 
expressions for P° and P°, involving terms up to first 
order in qP. Combining these results, we have 



(P^ 






In the position representation, P 



G{x) 



T- 



ip*{x)dxi^six) 
1 



(11) 
Oiq^)]Mx) dx. 

= —id/dx, we derive 

(12) 
1 



^{P:{k)+^P:{k))e-i'''' - ^P^Ak)+^P:{k)) 



where q = {f3/a)k. 

This expression can be combined with the previous re- 
sult, where we found |^s(x)p, producing the ratio 



i>*(x)dxipsix) _ dxi^sjx) _ G{x) 

tjj*{x)tljs{x) tps{x) li^sixW 



(13) 



which results in a differential equation for the wave 
function in the position space. In the general case in 
which ips{x) is a complex function, the knowledge of 
f{x) ~ G{x)/\ips{x)\'^ allows us to write a set of coupled 
differential equations 



dipljx) 

dx 

9V'™(a;) 

dx 



r{xWs{x)-r{xwr{^) 
r{xwr{x) + r{xWs{x), 



(14) 



where ipl'^'^{x) and /"^'""(a;) are the real and imaginary 
parts of ips{x) and f{x). The set of Eqs. (14) can be 
easily integrated with any numerical protocol given a 
known initial value, coming from the previous knowledge 
of iV's (2;) P, and the addition of an irrelevant global phase. 
On the other hand, following a similar procedure, the 
wavefunction in the momentum space can be measured 
by setting properly the control parameters to exchange 
the roles of X and P. Tuning the control parameters 
to produce the Hamiltonian H = Hg/SaxP, we can ob- 
tain an expression similar to Eq. ([5]) but in momentum 
representation, 17/1^(^)1^. To find ipsip)i we consider the 
Hamiltonian ^ but now with the condition (kX) <C 1, 
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FIG. 1. (Color online) Different steps needed to solve the wave function in position space of the state \ips) = (|0) + i|l))/-\/2 
using the method for pure states, a) Plots of P'^(k) + iP°{k) , which is real in our example (solid-black line), and the real 
and imaginary parts of P^(k) + iP°{k) (dotted-green and dashed-magenta lines). The inset represents the blow-up of the 
corresponding imaginary part, P°{k). b) Plots of |'i/)s(a;)p, (solid-black line), Re [G(a;)] (dotted-green line) and Im[G(a;)] 
(dashed-magenta line), c) Comparison between the real and imaginary parts of the wave function obtained using the method, 
(dotted blue and yellow lines) and the exact values calculated with the analytical expressions for {x\tps) = ({a^|0) +i{x\l))/\^, 
(solid-black lines), in this case the relation between a and /3 is /3/a = 0.5 x 10~^. The dashed green and magenta lines in (c) 
correspond to the results of the real and imaginary parts of tp calculated now with a less strict choice of P/a — 0.1, which 
shows the minor effects of the second order contributions of q in Eq (111 and, henceforth, the robustness of our method. 



and we find a differential Eq. for tps{p) analogous to 
Eq. (13). We present a pedagogical example in Fig. [I] 
In this manner, we have reconstructed the wave function 
of an arbitrary pure quantum state in the position rep- 
resentation by monitoring the evolution of an observable 
of the two-level probe. 



B. Mixed states 

In the laboratory, pure states can only be approximate 
and, in general, we deal with mixed states. In this case, 
a quantum system is described by a density matrix in- 
stead of a wavefunction. To have a deeper insight, in 
Fig. [2] we simulate numerically an experiment consider- 
ing the Hamiltonian H = hgacTxX. Note that in both 
examples, see Figs, [l] and [2J the dimensionless param- 
eters k and X are considered in a range where typi- 
cal experimental times are several orders of magnitude 
smaller than the involved decoherence ones [32 [33] . We 
compare then the wavefunction of a pure state with a 
similar one under the effects of a depolarizing channel 
£(|?/'s)(V's|) = (1 — e)|V's)(''/'s| +el/d. This channel maps 
the pure density matrix |?/'s)(i/'s| onto a superposition of 
this pure state and a fully mixed state I/d with proba- 
bility e, being d the dimension of the system. In Fig. [2] 
we show the robustness of the method for reconstruct- 
ing pure states in the presence of some decoherence pro- 
cesses. This numerical simulation was implemented by 
sorting a random number and comparing it with the the- 
oretical value of occupation probabiliy P°(fc), assuming 



1% of error in detectors. If the random number is less 
than the corresponding P|(fc), we assume the probe to be 
in the excited state, otherwise it would be in the ground 
state. Without the added error, we would approach the 
theoretical values for a large sampling. In the following, 
we extend our method to measure density matrices. 

Let us now consider the general case of the system 
described by the density matrix in the Fock basis Ps = 
'^Cn,m\'n'){m\. We will show that all matrix elements, 
Cn,m-, can be obtained using the same technique consid- 
ered in the pure case. The choice of /3 = in the one- 
dimensional version of the Hamiltonian M determines 
the evolution of the expectation value of the observable 
CT^, given by 



(a,)t-Tr[p,®|$)($|e^9*-^^a 



-igtcr^Xi 



(15) 



Performing both sets of measurements for the two initial 
probe states, we find 



P^ik) + t 



P°(^) = I]c„,„ [ dx 



^rmix)M^), (16) 



where k = 2gt and ipjix) — {x\j) is the wavefunction 
associated to the j-th Fock state. This last expres- 
sion yields a set of equations for the coefficients Cn,m- 
However, since the wavefunctions of the Fock states are 



real, Eq. ( 16 1 only provides partial knowledge about the 
density matrix. In consequence, we need to introduce 
more relations to complete a set of linearly indepen- 
dent equations allowing us to find all c„^m- To this end, 
we let the system evolve under the free-energy Hamil- 
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FIG. 2. (Color online) a) and b): Plots obtained using the described method for a state l-^s) ~ (|0) + |1) + |2) + |3))/2 in a 
depolarizing channel £{\tps){tp3\) = (1 — e)|V's){V's| +^I /d with e — 0.1. Dots show data from numerically simulated experiment. 
Green line is the fitting curve for the simulated data. Blue line is the expected result when e = 0. In c) the squared "wave 
function" obtained from data in a) and b) (green line) is compared with the expected |?/'s(a;)p for e = (blue dashed line). 



tonian H = hwa) a for a time t = Iq. This evolu- 
tion produces a change in the matrix elements of the 
system, i.e., Ps{to) = exp{—iuja^ato)psexp{iuja^ato) = 
Sn,mC„,mexp(-i(n- m)a;to)|n)(m|. Under this new 
initial probe state, Eq. (fl6|) becomes 



P^ik)+tP°{k) 






-i{n—rn)ujtQ 



X / dx e''-r„,{x)M^). (17) 



Choosing properly different values of the parameter tg, 
depending on the size of the density matrix p^, we obtain 
a set of linear equations allowing us the estimation of all 
matrix elements c„,„i. 

The proposed method for measuring the density ma- 
trix of a harmonic oscillator, via direct estimation of its 
matrix elements from Eq. (171, might be compared to 
some early Wigner function reconstructions |35l 136) . In 
those works, marginal distributions of the Wigner func- 
tion, associated to arbitrary field quadratures, are ob- 
tained using homodyne detection techniques. In our case, 
we encode the information of the density matrix elements 
in the level statistics of a probe that is coupled to the 
system with a Hamiltonian proportional to X. The ad- 
ditional phase rotation of the initial state, via the use of 
Hamiltonian H = huja^a, allows us to access information 
related to other quadratures, that is, linear superposi- 
tions of X and P. However, our method relies on probe 
measurements at arbitrary interaction times, making it 
difficult to trace any linear dependence with the involved 
Hamiltonians along the quantum evolution. 



III. PHYSICAL IMPLEMENTATION 

The methods described above could be implemented 
in different physical setups where two-level probes inter- 
act with quantum harmonic oscillators, as is the case of 



cavity QED and trapped ions in quantum optics, or cir- 
cuit QED in mesoscopic physics. The case of trapped 
ions can be easily implemented combining red and blue 
sideband excitations with suitable phases [57H39J . Here, 
we will give an example of how to construct the required 
family of Hamiltonians in the case of cavity QED. We 
consider then a field mode oscillator with energy huj and 
a two- level probe atom with transition energy Hujq. Let 
us assume that the system and the probe are in the 
joint state \'^) = |$)p (g) lips), where |$)p denotes the 
state of the probe and \4>s) denotes the state of the sys- 
tem to be measured. We allow them to interact via a 
Jaynes-Cummings coupling while the probe is driven by 
a phase-sensitive coherent field with frequency ujl- The 
total Hamiltonian describing this situation reads 00] 



H = ^a;o|e)(e| + hwa* a + hg{a* a + aa*) 



-hn{a 



^ —ii^Lt+if 



^iujLt—iip 



), 



(18) 



where a\a) is the creation (annihilation) operator of the 
oscillator mode and g is the coupling strength between 
the probe and the system. In the strong-driving regime, 
f2 3> {g, V'}: and under resonant conditions, uq = uj = 
lul, the Hamiltonian in the interaction picture reads 

i^int = y (|+^)(+^| - \-<p){-^\)iae-^^ + ah^^ 

= hga^X^ = hga^iaX + /3F), (19) 

where a^ — |+ip)(+ip| — \—<p){—ip\, hi the dressed basis 
|±^) = (Iff) ± e*'^|e))/y2, and X^ = {ae-'^P + a^e''f)l2. 
Therefore, we have shown how to build the desired Hamil- 
tonian, similar to Eq. (fl]). 



IV. SUMMARY 

In summary, we have developed a method for measur- 
ing the wavefunction of a quantum system in a pure state 



in position and momentum space, or the density matrix 
of a mixed state. This is achieved by suitable monitoring 
the evolution of a two-level probe coupled to the system 
via a Hamiltonian linear in position and momentum vari- 
ables. We expect that the proposed methods contribute 
to the already mature field of quantum state tomography. 
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